Berry-phase blockade in single-molecule magnets 
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We formulate the problem of electron transport through a single-molecule magnet (SMM) in the 
Coulomb blockade regime taking into account topological interference effects for the tunneling of the 
large spin of a SMM. The interference originates from spin Berry phases associated with different 
tunneling paths. We show that in the case of incoherent spin states it is essential to place the SMM 
between oppositely spin-polarized source and drain leads in order to detect the spin tunneling in the 
stationary current, which exhibits topological zeros as a function of the transverse magnetic field. 
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Single-molecule magnets (SMMs), such as Mni2 0,0 
and Fes [3|, |4| , have become the focus of intense research 
in the last decade since experiments on bulk samples 
demonstrated the quantum tunneling of a single mag- 
netic moment on a macroscopic scale. These molecules 
are characterized by a large total spin, a large magnetic 
anisotropy barrier, and anisotropy terms which allow the 
spin to tunnel through the barrier. Transport through 
SMM offers several unique features with potentially large 
impact in applications for magnetic devices based on 
SMM such as high-density magnetic storage as well as 
quantum computing applications Q. Recently experi- 
ments have pointed out the importance of the interfer- 
ence between spin tunneling paths in molecules and its 
effects on electron transport scenarios involving SMMs. 
For instance, measurements of the magnetization in bulk 
Fg samples (see Ref. |(|) have observed oscillations in the 
tunnel splitting AE S ,- S between states S z = s and — s as 
a function of a transverse magnetic field at temperatures 
between 0.05 K and 0.7 K. This effect can be explained by 
the interference between Berry phases acquired by spin 
tunneling paths of opposite windings using a coherent 
spin-state path integral approach 0, E| • 
The ability to manufacture molecular structures with 
fixed magnetic properties has given rise to a field known 
as molecular magnetism. To date, several experiments 
have shown the possibility to work with an individual 
SMM preserving the magnetic properties thereby 
demonstrating the Coulomb blockade effect at low tem- 
perature in a single SMM transistor geometr y 111! . The 
theoretically predicted Kondo effect in SMMs [l3,ll3| has 
not been observed yet. A theoretical description of the 
observed Coulomb blockade effect has recently been given 
by means of coherent spin states in Ref. • 
In this letter we propose the Berry-phase blockade effect 
by coupling an individual SMM to spin-polarized leads. 
We analyze the transport properties of the system in the 
Coulomb blockade regime for the ground state of a SMM 
in the presence of a longitudinal and transverse mag- 
netic field. Since the decoherence time between degen- 
erate spin states can be as low as T 2 — 10~ 9 s [lil llrij]. 



we work with incoherent spin states. We show that in 
the case of incoherent spin tunneling it is essential to use 
oppositely spin-polarized source and drain leads in order 
to be able to observe variations of the stationary current 
as a function of longitudinal or transverse magnetic field. 
In particular, the current can be suppressed due to the 
Berry-phase interference of the spin tunneling paths. In 
the case of fully polarized leads, complete current sup- 
pression coincides precisely with the topological zeros of 
the spin tunneling. Even in the case of partially polar- 
ized leads there are still fingerprints of the Berry-phase 
blockade which lead to significant changes in the station- 
ary current through the SMM. 

In the following we present our calculations that are valid 
for any SMM in the Coulomb blockade regime coupled to 
polarized leads and in the presence of a longitudinal and 
transverse magnetic field. We derive the (generalized) 
master equation for the low energy states and calculate 
the stationary current through the SMM for the cases of 
unpolarized, fully, and partially polarized leads. We ap- 
ply our results to the newly synthesized SMM N14, which 
has a spin of s = 4 and a ground state tunnel splitting of 
A£? S) _ S » 0.01 K at zero magnetic field correspond- 
ing to an angular frequency of uj = 10 9 s . 
Consider a SMM in the Coulomb blockade regime which 
is tunnel- coupled to two polarized leads at the chemical 
potentials fii, where I = L,R denote the left and right 
lead, respectively. The total Hamiltonian is given by 

H = H t +H s + H mi (1) 

where Hi = J2ika e ik c \ka c ik<? represents the energy of the 
leads. c lk creates an electron in lead I with orbital state 
k, spin cr, and energy eik- The coupling of the leads to 
the molecule is described by H m = J2i P ka t( p c ikadp<T + 
H.c., where tf p denotes the tunneling amplitude and d pt7 
creates an electron on the molecule in orbital state p. The 
term 7i sp m is typically given by the spin Hamiltonian of 
a SMM in an external transverse magnetic field H± and 
a longitudinal magnetic field H z , i.e. 
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where the easy axis is taken along z, S qt ± 
and the integer index q denotes the charging state of the 
SMM, where q = — 1 adds one electron to the molecule 
and q = when the molecule is neutral. The transverse 



magnetic field Hj_ = H x + iH y = \H±\e t(p lies in the 
xy plane. We use the abbreviation h±_ = gpsH±. In 
this Hamiltonian, the dominant longitudinal anisotropy 
term creates a ladder structure in the molecule spectrum 
where the eigenstates | ± m q ) of S z are degenerate. The 
weaker transverse anisotropy terms couple these states. 
The coupling parameters depend on the charging state 
of the molecule. For example, it is known that Mni2 
changes its easy-axis anisotropy constant (and its total 
spin) from A„ = 56 peV (S = 10) to A_i = 43 peY 
(S-i = 19/2) and A_ 2 = 32 ^eV (5__ 2 = 10) when singly 
and doubly charged, respectively 1 1 Qfj . Experiments with 
NU show that B4 i9= o = —0.003 K, i.e. is negative 



|18| . In this case, in order to see the Berry phase oscilla- 
tion, a magnetic field H± must be applied in the xy plane 
[lo| along specific angles ip(B q ,B^, q ). It is also possible 
to tune the tunnel splitting by means of the longitudinal 
magnetic field H z . 

For weak coupling between the leads and the SMM we 
use the standard formalism suitable to describe a system 
(SMM) coupled to a reservoir (polarized leads) 20] . The 
master equation describing the electronic spin states of 
the SMM is given in Born and Markov approximation by 



oscillations. In this case unpolarized leads can be used 
to measure the tunnel splitting between the coherent spin 
states (|s') + |-s'))/V2 and (|s') - \-s'))/V2 by varying 
the gate or bias voltage. However, only partially or fully 
polarized leads allow us to probe the incoherent tunnel- 
ing rate r Sj _ s between the ground states s and — s for 
q = and also between s' and — s' for q = —1. As we 
prove below, both incoherent tunneling rates r s _ s and 
Ts/^s' contribute to the total polarized current through 
the SMM. The sequential tunneling rates for absorption 
of an electron in Eq. (j3J for ground states with spin ,s 
and s' and energy differences A± s ^± s in the case of low 
temperatures are given by 



(4) 

and the tunneling rates for the emission of an electron 
are given by 
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w- a ,-* = ZiWLl_j,w® i _, = w®[i-MA- a ,-.,)] : 

(5) 

where fi(A s , jS ) = [l+ e ( A >',>-^/ kT ]- 1 is the Fermi func- 
tion. The diagonal elements of J3J| yield 



p s = -[p, H] s , s + ^2 PnWs,n ~ Ps^2 ^™> s > 
and the off-diagonal elements of Q yield 



(6) 
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w mJ l m, 71 P 771,11 ) (3) 



where 7 m ,„ = \ J2i(Wi,n + Wi, m ) + 1/T 2 is the total de- 
coherence rate which contains the spin decoherence time 
T 2 due to e.g. nuclear spins and the rates W m>n of tran- 
sition between the states of the SMM. Figs. (JTJ and 
show the ground states s and — s of a SMM placed be- 
tween unpolarized and polarized leads, respectively, w^j 
represents the spin-dependent transition rate from the 
I = L,R lead to the SMM and are defined in Fermi's 
golden rule approximation by w^p — 2ttDv^ \t^\ 2 /h and 

= 2irDv^ \t^ \ 2 /h, respectively, where D is the den- 
sity of states and in and in' are fractions of the num- 
ber of spins polarized up and down of lead I such that 



,(0 



1. 4° and t) l > are the tunneling amplitudes 



T 



of lead I, respectively. Typical values for the tunneling 
rate of the electron ran ge f rom around w = 10 6 s" 1 to 
w = 10 10 s" 1 (see Refs. £l|2l|). 

In order to see coherent spin tunneling, the de- 
coherence time must be increased for example by 
deliberation such that I/T2 -C AE s ,- s /h and at the 
same time the contact to the leads must be so weak that 
W m . n <C AE Sl - s /h. Another possibility is to increase the 
transverse magnetic field \H±\ beyond the Berry-phase 
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Since we are interested in the long time behavior t 3> 
l/7m,n, we can set p SlS i — in equation J7J| to obtain the 
following coupled differential equations for the diagonal 
elements of the density matrix 



Ps 



AE, _. 
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2h J {H s -H- s ) 2 /fi?+i 2 g _ a (P- s P s } 

+W s , s >p s > - W s , s p s , 



P-s = 



( AE S - B \ 2 2 7s _ s / _ \ 

\ 2h J (H_ 3 -H 3 ) 2 /h 2 +~fl_ a P~ s > 



+ W- S -s'P-s' - W- S '- S p- 



(8) 



(9) 



The other two differential equations are obtained by just 
replacing s «-> s' in the above equations. Solving the set 
of differential equations for p s , p s > and in the 
stationary case (t ^S> 1/W mi „) we obtain 

Ps = {w StS .{W- s ,,- s + r Sj _ 5 )rv,_ s , 
+w_ s ,_ s /(w s / )S + r s ,,_y)r S; _ s )/7 ? , 

Ps = (w a . s ,{W-s.-s> +r fl ,,_ s or s ,_ fl 
+w_ ai _./(w,/,. + r, i _ i )r y ,_^)/t J) 

Ps' - {W s ,,s{W-s.-si + T s ..-s-)T s ,-s 

+w-s>,- s (w s y + r s ,_ s )r s ,,_ s o/r7, 
Ps- = (wv, a (w_ a / ) _ s + r a ,_ a )r a /,_y 
+w-s>,- s (w s y + r a /,_ a /)r a ,_ a )/r ? , 



(10) 
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where r\ is a normalization factor given by 



V = (T4,-4 + r 7/2 ,-7/2)(l / l / -4,-7/2W / 7 /2,4 

+1^4,7/2^-7/2,-4) + 2[r 4 ,-4^_ 4 

+r7/2,-7/2M / -7/2, -4^7/2,4 + r4,-4r 7 / 2 ,-7/2 
(WV/2,4 + V^_7/2,-4 + ^4,7/2 + W_ 4 _ 7 / 2 )], 



7/2^4,7/2 



(11) 



such that J2 n Pn 
{AE S 



1. The incoherent tunneling rate is 



V 2h 



(12) 



which is a Lorentzian as a function of the longitudinal 
Zeeman splitting H s — H- s = gfJ,sH z [s — (— s)]. We now 
proceed to define the current through the SMM in terms 
of the density matrix elements. In the case of Ni 4 we 
have s = 4 and s' = 7/2, therefore the current reads 



e (W / 4 ! 7/ 2 P7/2 + ^-4,-7/2^-7/2) 



(13) 



Taking into consideration the asymmetry of the leads ,i.e. 

„J L ) -1 „„(«) 



at 



7^ w\* , and restricting ourselves to the case of un- 



polarized leads, i.e. 



,{R) 



T 



= 1/2, we 



obtain the following conditions for the transition rates 
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W_ 



7/2, -4> 



it; 



4,7/2 



4,-7/2- 



(14) 



Substituting the values of p 7 / 2 and p_ 7 /2 into Eq. 113|l 
and using Eq. I|14|l we obtain 



(15) 



■lip 



ir- 



7/2,4 



4,-7/2 



which does not depend on the tunnel splitting energy of 
the SMM. Thus it is impossible to observe Berry-phase 
oscillations for the case of unpolarized leads and inco- 
herent spin states. Eq. (|15fl can be interpreted as two 
resistances in series |17J, where the only transitions that 
contribute to the current through the SMM are 4 <-> 7/2 
and -4 <-> -7/2 (see Fig.HJ. 

In the case of leads that are fully polarized in opposite 



directions, i.e. v 



(L) 



,(R) 



1 or v 



(L) 



(R) 1 
v; ' = 1, we 



get the following conditions for the transition rates: 



-4,-7/2 = ^7/2,4 = 0, W<4,7/2 = W-j/2 ,-4 = 0. (16) 



Choosing the case 



(L) 



,{R) 



1 and using Eq. l|ltj|) 



we can then substitute the values of p 7 /2 and P-7/2 into 
Eq. (O to obtain 



e 

T, 



w. 



7/2,-4 



r 



4,-4 



4,7/2 



r 



7/2,-7/2 



(17) 



which reflects the fact that the current through the SMM 
depends on the tunnel splittings and can be interpreted 
as four resistances coupled in series in a loop (see Fig.|5J). 
Notice the clockwise direction of the transition rates be- 
tween the different states s, s'—s' and — s of the SMM. 




FIG. 1: Diagram for the transitions between the ground states 
4 «-» 7/2 and —4 «-» —7/2 in the case of unpolarized leads. The 
transitions arise from the sequential tunneling of unpolarized 
electrons in and out of the SMM. E c is the charging energy. 




FIG. 2: Diagram for the transitions between the ground states 
4 — > 7/2 — ► —7/2 — > —4 in the case of fully polarized leads in 



opposite directions v\ 



(R) 

I 



If we chose to work with fully polarized leads of the form 



,(R) _ 



1 , then the direction of the transition rates 



between states would be opposite, i.e. anticlockwise. 
Fig. (PJ shows the current as a function of the trans- 
verse magnetic field H± for fully polarized leads. If the 
tunnel splitting Ai? 4i _4 or A_E 7 / 2 ,-7/2 is topologically 
quenched, then r4 : _ 4 or r 7 / 2 .-7/2 vanishes [see Eq. l|T2)l]. 
which leads to complete current suppression according to 
Eq. (|1 T|) . Since this current blockade is a consequence 
of the topologically quenched tunnel splitting, we call it 
Berry-phase blockade. Note that the current can also be 
suppressed by applying the longitudinal magnetic field 
H z which follow immediately from Eqs. (|12f> and 117(1 . 
If we consider now partially polarized leads (i.e. vX 



(R) 



> 



,(«) „( L ) 



■'I 



> v 



(L) 



) and calculate the current through the 



4 





FIG. 3: The graph shows the logarithm base 10 of the sta- 
tionary current versus the transverse magnetic field H± for 
v\ R) = v{ L) = 1, v{ R) = u\ L) = 0, w{ L) = 1 x 10 9 *- 1 and 
w l = 10™j L '.The polarizations of the left and right lead are 
given by P R = -P L = vf - vf = 100 %. At the zeros of the 
tunnel splitting AE a ,- a or A-E s ',-s' the current is completely 
suppressed. The scale varies from I p = 0.1 nA to 1 fA. 



FIG. 4: The graph shows the logarithm of the stationary 
current versus the transverse magnetic field H± for vi R ^ = 

v{ L) = 0.8, v{ R) = v\ L) = 0.2, w[ L) = 1 x lO 9 ^ 1 , and w{ R) = 
Ww, . The polarizations of the left and right lead are given 
by P R = -P L = i/f - i/f = 60%. In this case at the zeros 
of the tunnel splitting AEs,~s and AE s > t - s i the current is 
suppressed approximately by a factor of 3. 



SMM, we obtain 

Ipp = e[r 4 ^4VF_4 i _ 7/2 VF4 i7/2 (W / _ 7/ 4 ! _4 + W V 2 A ) 
+r 4j _4r 7 / 2 , -7/2(^-4,-7/2 + W / 4 !7 / 2 )(W / _ 7 /4,_4 + Wr/2,4) 
+r 7 /2.-7/2W A _ 7 /4.-4W 7 7/2,4(W / -4 -7/2 + W 4 j/ 2 )}/V- 

(18) 

Fig. 01 shows the current as a function of the transverse 
magnetic field H±_ for partially polarized leads. The 
Berry-phase blockade is still visible in the stationary cur- 
rent even at a spin polarization of 60 %, which makes it 
experimentally accessible since recent experiments have 
achieved near 100 % spin polarization |23| . 

In summary, we have shown the Berry-phase blockade for 
a SMM placed between polarized leads. This behavior is 
due to Berry-phase interference of the SMM spin between 
different tunneling paths. We have shown that in the case 
of incoherent spin states it is essential to use polarized 
leads in order to observe the Berry-phase blockade. 
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